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Abstract 

Pure leptonic and semileptonic rare B decays, Bs ^J-^ and B — > Xsfi^ fJ.~ in the minimal super- 
symmetric Standard Model (MSSM), in particular, gluino and neutralino contributions to the decays, are 
discussed under Br(_B —+ Xs^) and other experimental constraints . The general scalar quark mass matrices 
as the new sources of flavor violation are considered. We present Wilson coeflicients of bs transitions from 7, 
Z, and neutral Higgs boson penguin diagrams by using vertex mixing method to deal with scalar down-type 
quark flavor changing and also give their expressions in MIA to show different sources of enhancements. 
We find that under the experimental constraints, with large mixing of left-handed and right-handed sbot- 
tom, C[q can be enhanced by 10% compared with SM, in two cases, heavy gluino and fine-tuning between 
523 ^ {623'^) and -^(^If^) terms in cj''. Particularly, Cio and C[o can reach a 20% enhancement in some 
regions of parameters under experimental constraints. When CP-odd Higgs A" is not too heavy (~ 250GeV), 
and tan/3 is large (~ 40), neutral Higgs boson penguins with gluino and down-type squark in the loop can 
significantly contribute to the bs transition and the contributions can compete with those due to the chargino 
and up- type squark loop. 
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1 Introduction 



Pure leptonic and semileptonic flavor changing neutral current (FCNC) rare B decays, 

^ M+M" , B ^ (1) 

have received much attention recent years due to clear backgrounds and ongoing experiments at BaBar [1] 

and BELLE and forthcoming projects at Tevatron [2] and LHC [3] as well as sensitivity to models beyond the 
Standard Model (SM). The current experimental result of Br(i? — > XgfJ.'^fJ.^) by the BELLE collaboration [4] is 

Bv{B X,M+/x-)exp = 7.9 ± 2.lt\\ X 10-^ (2) 

and CDF [5] upper limit on Br(Bs is 

Br(Bs ^ M"^M")exp < 2.6 X 10"^ at 90%C.L. (3) 

In the SM, these processes vanish at tree level, while they occur at one-loop level with the charged gauge 
boson and up- type quarks in the loop. In the minimal supersymmetric Standard Model (MSSM), there are 
five kinds of contributions to partonic level process b s^'^ at one-loop level, depending on specific particles 
propagated in the loop, (1) Standard Model gauge boson and up-type quarks (SM contribution); (2) charged 
Higgs iJ^ and up- type quarks (charged Higgs contribution); (3) chargino and scalar up-type quarks (chargino 
contribution); (4) neutralino and scalar down- type quarks (neutralino contribution); (5) gluino and scalar down- 
type quarks (gluino contribution). The flavor structure of the sfermion sector in MSSM depends on the soft terms 
which are determined by the supersymmetry breaking mechanism, in addition to the superpotential. In the 
minimal flavor violation (MFV) scenarios of MSSM, squarks are assumed to rotate in flavor bases aligned with 
the corresponding quark sector and the only source of flavor violation is the usual Cabibbo-Kobayashi-Maskawa 
(CKM) matrix in SM. In some MFV scenarios such as the constrained MSSM (mSUGRA, string-inspired flipped 
SU(5), etc.) and gauge mediation supersymmetry breaking (GMSB) where the soft terms at some high scale 
(the grand unification scale or Plank scale or messenger scale) are characterized by the universality of sfermion 
masses and the proportionality of the trilinear terms, the flavor violation in sfermion sector at the electroweak 
(EW) scale is generated radiatively and consequently, in general, small ^. Therefore, comparing with the first 
three kinds of contributions, the last two kinds of contributions, i. e., neutralino and gluino contributions, are 
negligible [13]. 

There are new sources of flavor violation in MSSM. Besides the CKM matrix, the the 6x6 squark mass 
matrices are generally not diagonal in flavor (generation) indices in the super-CKM basis in which superfields 
are rotated in such a way that the mass matrices of the quark field components of the superfields are diagonal. 
This rotation non-alignment in the quark and squark sectors can induce large fiavor off-diagonal couplings 
such as the coupling of gluino to the quark and squark which belong to different generations. There exist two 
different kinds of methods to deal with flavor changing vertices induced by fiavor mixing in the squark mass 
matrices in the literature [9]. One works in quark and squark mass eigenstates with induced fiavor changing 
couplings, so called "vertex mixing" (VM). The other method, "mass insertion approximation" (MIA) [10], 
works in flavor diagonal gaugino couplings gqq and diagonal quark mass matrices with all the fiavor changes 
rested on the off-diagonal sfermion propagators. The MIA can be obtained in VM through Taylor expansion of 
nearly degenerate squark masses rrig. around the common squark mass rriq, m? ~ m,?(l + A^). Thus MIA can 
work well for nearly degenerate squark masses and, in general, its reliability can be checked only a posteriori. 
However, for its simplicity, it has been widely used as a model independent analysis to find the constraints on 
different off-diagonal parts of squark mass matrices from experiments [11]. It is clear that VM remains valid 
even when fiavor off-diagonal squark mass matrix elements are large and there is no approximation which has 
been assumed. 

Because the experiment of Br(B — > Xg'y) only constrains \CT{mh)\^ -\- \C'rj{mh)\^, with the overall sign flip of 
Crifnt) from MSSM vs SM, the branch ratio of i? ^ Xsl~^l~ can be enhanced from the term Re(C7(mh)Cg^'''*(m(,)) 
in the constrained MSSM with gluino and neutralino contributions neglected [12, 13, 15], while in these analysis 

^In some large tan/3 regions of the parameter space it becomes non- negligible [14]. 
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supersymmetry contributions to Cg and Cio are at most changed by ±5% compared with the SM values. The 
chargino contribution in the extended MFV model is analyzed in ref [16]. B — > Xsl~^l~ has been analyzed in the 
left-right supcrsymmctric model recently [17]. In MSSM, gluino induced FCNC process b is studied 

in ref [18], and chargino, gluino induced effects are studied in ref [19] in MIA. In these works the contributions 
from exchanging neutral Higgs bosons (NHBs) are not included. However, for I = /i,r, when tan/3 is large 
and is not too large (say, 250 GcV), the NHBs contributions can become significant due to the tan'^/3 
enhancement of the corresponding Wilson coefficients in some regions of the parameter space [20, 21] and the 
NHBs contributions to B — > XJ~^l~ and Bg — > [q = d,s, I = fi, r) in the constrained MSSM have been 
investigated in refs. [20, 21, 22, 23]. Using the VM method, the gluino and ncutralino induced FCNC processes 
b sjX^ iJL~ [24] and Bq jJ''^ IJ'~ [25] have recently been analyzed, including the NHBs contributions, in 
the constrained MSSM and in MFV models respectively. In SUSY models with non-minimal sources of flavor 
mixing, the constraints on different flavor violation parameters from Br(_B Xg^) have been considered in 
ref. [26, 27, 28], and Bs,d at large tan/3 has been investigated in ref. [29]. 

It is well-known that the effects of the primed counterparts of usual operators are suppressed by ^ and 
consequently negligible in SM because they have the opposite chiralities. In MFV models their effects are also 
negligible, as shown in ref. [25]. However, in MSSM their effects can be significant, since the flavor non-diagonal 
squark mass matrix elements are free parameters. Part of the primed counterparts of usual operators relevant 
to B rare leptonic and semileptonic decays have been considered in ref. [19]. 

In this paper, we extend our previous analyses to include gluino and neutralino contributions and all opera- 
tors responsible for B rare leptonic and semileptonic decays in MSSM. We calculate the Wilson coefficients using 
the VM method and also give their expressions in MIA to show different sources of enhancements. In numerical 
analyses we take into account constraints from Br(i3 ^sl), ^Mb^ and the lower bounds of superpartner 
masses and Higgs masses as well as i? ^ X^g and hadronic charmless B decays. We have carefully analyzed dif- 
ferent sources of enhancements of c|q (arising from the Z penguin) and Cq| ^ (arising from the NHB penguins), 
related to the general scalar-down quark mass matrix. We find that under the experimental constraints, with 
large mixing of left-handed and right-handed sbottom, C^q can be enhanced by about 10% compared to the 
Standard Model in two cases, heavy gluino and fine-tuning between ^23^(^23^^) and ^23 ^(^23^^) terms in C^'^ . 
In particular, Cio and C'iq can reach a 20% enhancement in some regions of parameters under experimental 
constraints. When CP-odd Higgs is not too heavy (~ 250GeV), and tan/3 is large (~ 40), neutral Higgs 
boson penguins with gluino and down-type squark in the loop can significantly contribute to the bs transition 
and the contributions can compete with those due to the chargino and scalar up-type quark loop. 

The paper is arranged as following. In section 2, we define our notations and consider the effective Hamil- 
tonian and branching ratios of pure leptonic and semileptonic rare B decays. In section 3 we briefly rcicall the 
squark mass matrices and discuss the choice of parameters. In section 4, we present our numerical analysis on 
the possible enhancement of C'['q and Cg| ^ in the case of switching on only the gluino (or neutralino) and SM 
contributions. We search for maximums of Cio and C[q under experimental constraints, switching on all the 
contributions. Section 5 is devoted to give the numerical results for B Xsl^l~ and Bg Z+Z~. Conclusions 
are drawn in section 6. In the Appendix, Wilson coefficients at scale are given. 



2 Effective Hamiltonian 



The effective Hamiltonian for B XJ^l and Bg can be written as 

" 6 10 10 

5]a(M)0,(M) + 5^(Q(M)0,(/x) +CKm)OUm)) +E(CQ.(M)Qi(M) + t?Q,(M)QUM)) 



V2 



i=7 



(4) 



where At = VtbVf*^, Oi and Qi (i=l,...,10) can be found in ref. [30] and [31] respectively, and the primed opera- 
tors, the counterpart of the unprimed operators, are obtained by replacing the chiralities in the corresponding 
unprimed operators with opposite ones. The explicit expressions of the operators governing B XJ~^l~ and 
Bs ^ are given by 



O7 



167r2 



mb{sa^.PRb)F^-', O'r 



167r2 



mb{sa^,uPLb)F'' 
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2 2 
2 2 

Qi = ^{sPRbm, Q'i = ^{sPLb){ll) 

2 2 

Q2 = ^{sPRb){Hl), Q'2 = ^{sPLb){hnl) (5) 

where = (1 — 75)/2 and Pr = (1 + 75)/2^. We also consider the operators 

Os = J^m6(s„T«^<7^,PK^);3)G°''^ = Jl2™b(s„T«^a^,Pi6^)G»'''' (6) 

in order to include the constraints from B — > X^g and hadronic charmless B decays into the analysis. In SM the 
Wilson coefficients of the prime operators are suppressed by ^ with respect to those of unprimed operators and 
the statement is also true in MFV scenarios of MSSM [25]. However, in MSSM, the statement is, in general, not 
valid due to the presence of new sources of flavor violation. The running of Wilson coeflicients Cj and CQ^ from 
m-uj to mt, in the leading order approximation (LO) is given in rcfs [30] and [31] respectively. The evolution of 
part of the primed operators has been given in ref. [26] . Although the mixing between Oj in the next- to- leading 
order (NLO) has been studied, the mixing of Oi with Qi in NLO has not been given. So we shall use only the 
LO results for consistence. We present raw scale Wilson coefficients and Cq. in Appendix. In order to see 
the dependences of Wilson coeflicients on new flavor violation parameters, as an illustration, we present all the 
relevant Wilson coeflicients at the scale induced by gluino in MIA as follows 

^ ^ ^ rUg ■^'^ rUg 



+ 



Cf\m,) = -0.21(i^^)24"(«^)+0.17(^^^)^4^«(«^)<53^3^« + 

0.675i^''^''''\6ity - 1.685i^''^''''^6i^'' (9) 
C'fo'\mb) = -8.374^^(''''^(5#^)2 + 21.04''''(''''^5#^ (10) 

Cq^ (m,) - 4.61( ) ( 40 ) 400GeV'^23 

cS:\mt) = T<^(m,) (12) 

Where = ^"^'"^^^'C,'^'^'-''^^^ (M|^b is the scalar quark mass matrix, q = u,d, A,B = L,R, i,j = 1, 2, 3, 
m is the common scalar quark mass), the one- loop functions (which are given in appendix) at a; = — | = 1 have 
been used, tan/J is the ratio of two Higgs vacuum expectation value vu and vd, tan/3 = and niHO = ^a", 
which is a good approximation in the case of m\o/'m% ^ 1, has been assumed. Here we have expanded the 
Wilson coeflicients to the double MI, as investigated in ref [27] , which is non-neglectable if the mixing between 
left-handed and right-handed sbottoms is large. In eqs.(7) and (8) the flrst two terms are suppressed by ^ (as 

^Note that Qi,2 = ^ ^^^^ Os,p and Os^p are used in some papers (see, e.g., ref.[32]). 
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can be seen from Appendix, although they are enhanced by a factor — , they suffer from the suppression). 
In eq.(9) for Cg (mi,), the first two terms come from the 7 penguin, which is suppressed by But eq.(lO) 

g 

for cfJ(TO6) and the last two terms in eq.(9) from the Z penguin is not suppressed, which is noticed as non- 
decouphng of the bsZ couphng in ref. [33]. Therefore, there is a possibihty that Cio can be enhanced even 
if the gluino is heavy. In eq.(ll) and eq.(12) for Cq^ ^, we present 62^^^^^^ contribution in the MIA, which 
is enhanced by a factor Compared to the chargino contribution, the gluino contribution is the same 

important provided that mi, tan /3 ^ mg. We can read that Cq| ^ can be large if CP-odd Higgs A'^ is not too 
heavy (say, < 400 GeV) and tan /? is large. 

The leading order Bs — > Xs'y branching ratio normalized to Bt{B Xcev) is given as 

Br(B ^ X,7) = ^l^l'Br(i? - X,ei>)(\C^{mi,)\'' + \C',{m,)\^), (13) 

where \/z = rn^°^^ / m^°^^ , f{z) is the phase space function. 

The branching ratio (Br) B — > normalized to Br(B — > XcCu) is given as 



dBr(5^X.,V-) ^ BriB^X^e,)^il-srJl-'-^^^Dis), 



ds ' ''i-K^f{zy ' V s 



cb 



D{s) = \Cf {mt)\\l + —){! + 2s) + A\Cr{mt,)f{l + —){! + -) + 

s s s 

|Cio(m6)n(l + 2s) + —(1 - 4s)] + 12Re(C7(mb)Cf *(m6))(l + — ) + 
s s 

l\CQ,{mh)\'{s~4f) + ^\CQ,{mb)\^s + 6MCio{rnb)C*Q^{mb))t + 

iCiimt,) ^ C'iinn,)). (14) 

where i = ^ and s = — ^jr— • In SM, Crimb) and Cf^*{mb) have opposite sign, so the term Re{C7{mb)Cf^* (mb))) 

decreases the overall results above. In MSSM, there exist some parameter regions, where supcrsymmetry con- 
tributions to C7 can have opposite sign compared with CY^{mb), and even make the sign of Ci{mb) opposite to 
but its size approximately equal to Cf^'(m{,). In this case the term Re(C7(mb)Cg^*(m6))) adds construc- 
tively to dBr{B^X^^+^.-) ^ ^i^j^j^ 

can make the Br of B — > X^Z+Z enhanced by approximately 50% compared 

to that in SM [15]. Prom the formula above, we can see that with large Cq| ^, S — > XsH~^fj,~ can be enhanced 
greatly. 

The branching ratio Bs — > n'^n~is given as 

Br(B,^/zV-) = ^f^^^LrsjIJAtlVl - 4m2[(i _ 4m2)|CQ,(m,) - C'Q,{mb)\^ + 

\CQ,{mb) - C'Q^imb) + 2m(Cio(m6) - C[o{mb))\^] (15) 

where m = m^/niBs- With large Cq| ^, Br(Bs — > can be enhanced by several order of magnitude. 



3 Squark mass matrices 

The 6x6 squark mass-squared matrices in the super-CKM basis have the structure 



M- M- 

= I mT Mr ) (16) 



^The (Sj^^'^^^ contributions to Cq^ from self-energy type and penguin diagrams are cancelled provided that m^o = m^o 



^^.ivii>^u^.^.io ^Ql,2 

and MIA is valid 
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where q = u(d) represent the up (down)-type squark. Note that differently from Mi, ^ is not 

hermitian. In general, A^|^^ ^„ are off-diagonal. The 3x3 submatrices are given by 



i^LJ^j = KJy + + cos2/3mi(/^ - sin" (?v^)) 

(■^Lh)*^ = (™L)»i + (™L + '5? cos 2/3m| sin^ 0w) 

iMl^^)ij = vuAuij - SijUiuiH cot P 

(-^dL^v = voAdij - SijiTLdiiJ.ta.TLP (17) 

where niq^ with g = d, w is the quark mass of generation i, and Qg are the third component of weak isospin 
and electric charge of quark q respectively, /x is the Higgs superfield mixing parameter, Auij, A^ij are trilinear 
higgs-squark-squark coupling. 

Because of SU(2) gauge invariance, the squark mass matrix to|^ is intimately connected to via 



l^CKM"^iVbKM (18) 



Then, from m\ , we can get m? , and vice versa. 

dL 

Furthermore, we assume for the sake of simplicity that there are no new CP-violating phases, besides the 
single CKM phase. Thus, we in general have twenty seven new flavor violation parameters totally from squark 
mass matrice. However, only five of them are involved in the transition b to s in our analysis (see below). 

Because we concentrate on the b — s transition, only left-left, right-right and left-right 2 — 3 mixing terms 
are directly dependent. In order to simplify the analysis we shall keep only these 2 — 3 mixing terms non-zero 
and set all the 1 — 2 and 1 — 3 mixing terms to 0. We also keep the third generation 3 — 3 left-right mixing 
term non-zero. The first generation 1 — 1 and second generation 2 — 2 left-right mixing terms are set to for 
simplicity. We parametrize the non-vanishing 2 — 3 off-diagonal term as 

{Ml^j2s = S^i^^iMlM^lLLh^ (19) 
In the super-CKM basis the fields qu, qm {i = 1, 2, 3) are related to the mass eigenstates (a = 1, 6) by 

qL,R = rt^ ^jg, (20) 

Where the matrix F^^ ^ (g = {/,£>) is a 3 x 6 mixing matrix. 

In order to simplify the discussion further, we assume all the diagonal elements of squark mass matrices are 

equal to a commen SUSY scale at electroweak scale, niq — rriu^.. = rriuj^.. = wj^., = J^ij^,, (i = 1,2,3). Then 
we have three flavor conserved parameter rriq, S]^^^, ^33^ a-^d five flavor violation parameters involved in the b 
to s transition, 623^, ^23^, 623^, ^23^^, (^23^^) from the squark mass matrices. In addition to these parameters, 
we have also the following free parameters: gaugino masses Mj (i=l,2,3), CP-odd Higgs boson mass M^o, fj,, 
and tan/3. 



4 Numerical analysis of Wilson coefficients Ciq and Cq\ ^ 

As analyzed in section 2, C[q from the Z penguin and Cg| ^ from neutral Higgs boson penguins can be signifi- 
cantly enhanced. In this section, we present our numerical analysis of Wilson coefficients C'[q and Cq| ^ versus 
fiavor violating parameters under experimental bounds, particularly Br(_B Xs"f), in the case of switching on 
only the gluino (or neutralino) and SM contributions in the first three subsections (fourth section), and then 
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switching on all the contributions from the W boson, charged Higgs, chargino, neutralino, and gluino in the 
last subsection. We show there are two cases, the heavy gluino, and fine-tuning (i. e., to finely tune 62^^^^^^ 

and $2^^^^^^ terms in C^'^ makes the constraint from b ^ satisfied) to enhance C[q. In general Cg\'mb) 
in MSSM is enhanced by at most five persent compared to SM. We can see the reason why it can not have a 
large enhancement from the expressions of Cg'^ in the Appendix. The gamma penguin and box contributions 
to C^'' are suppressed by a factor of m'^/m'i± -0 -■ And the Z- penguin contributions to Cg'^ are suppressed by 
(— 1 + 4,s^,,). So we shall not discuss its dependence on new flavor violation parameters in most of part of the 
numerical analysis hearaftcr. In our numerical analysis we use the expressions of Wilson coefficients obtained 
by the VM method, i. e., those given in Appendix. In order to simpliiy our notation and express new physics 
effects we introduce the following quantities, 

where C'W'"^^™ is the Wilson coefficient of the operator O^P_ in MSSM and Cf " in SM, and i = 7, 8, 9, 10. 
The experimental measurements of mass differences in — system give the following bound [34] 

AMs > 14.4ps-^ . (22) 

We consider the constraints of AA'/s on (5's as in ref. [SG]**. Current c;xperimental measurements of inclusive 
decay B X^^i at ALEPH [6], CLEO [7] and BELLE [8] produce the world average value 

Br(B ^ Xs7)exp = (3.23 ± 0.41) x lO"'' . (23) 

We use 2cr Br(B ^s7) bound in our numerical analysis. 

4.1 Heavy gluino contributions to C^l 

As we have noticed in section 2, though C7 ' has ^ enhancement, it still suffers from ^ suppression. Then 

C^'"* is suppressed by as presented in Eq.(7). Nevertheless the Z penguin contribution to c{q is non- 
decoupled when nig is large. Therefore, we expect that when the gluino is heavy, the b sj constraint can be 
easily satisfied and a large enhancement of C[2 can occur. Including, in addition to the SM contribution, only 

the gluino contribution, we show 623^, ^23'^; ^23^^; and 623'^ dependences of c[2 respectively, i. e., in each 
case only one non-zero off-diagonal parameter enters, in Fig. 1, where Mq = Mj = 800GeV, Mi = lOOGeV, 
M2 = 1200GeV, M3 = 3000GeV, 11 = 3200GeV and tan/3 = 50. The choice of large n and tan/3 is to ensure 
large left-handed and right-handed sbottom mixing, ^33'^ = —0.75. 

In Fig. l.a, we give i?io as a function of 623^- When 62^^ is near 0.6, i?io is the order of 10%, where the 
sign of C'j{mb) is flipped, and the lightest sbottom mass is near the low bound obtained in experiments. We 
show i?io as a function of 82^^ in Fig. l.b. Rio can also reach the order of 10%. 

When we switch on only S!^^^ or 523^^, we assume some mechanism to render C-j{m}j) to be 0, e.g., some 
non-zero (^23^ or (^23^ contribution to Cr{myj) nearly cancells the SM contribution. We present R'^q as a 
function of ^23^^ iii Fig. I.e. When ^23^^ is near —0.6, R[q can be as large as 7%. R'iq as a function of S23'" are 
shown in Fig. l.d. When 623^ is near —0.15, R[q can be as large as 7%. 

4.2 Fine-tuning of ^23^^^^^ and 62^^^^^^ terms in C^'^ and gluino contributions to 

^10 

It is obvious from Eq.(7) that C^'' can be finely tuned to zero with the large cancellation between ^23^^^^^ and 

623^^^^^ terms. However, in this case the ^23^^^^^ and ^23^^^^^ contributions to C['q can be large. 

■^Thc analysis in ref. [36] is based on the old experimental bound, AMg > 15.0ps~^ [35]. An update analysis in MSSM is in 
progress [37]. 
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We show the correlations of 623^ vs 623^ and Rg vs Rio in Fig. 2. a, where fine-tuning of 623^ and 623^^ 
terms in C7 have been carried out. The other parameters are Mq = SOOGeV, M3 = SOOGeV, /x = 1200GeV, 

and tan/3 = 50. Large /i and tan/? is to ensure large sbottom mixing 

<^33 ^ = -0.72. Wc find that the largest 
Rio is 40%. However, when we impose the constraints from B X^g and hadronic charmless B decays which 
require |i?8| cannot be larger than 10 [38, 16], Riq can only be as large as 8%. 

In order to analyze S23^^^'^\ wc set CT{mi,) to and keep Cg(TOh) and Cio{m(,) to be their corresponding 
SM values. We make the fine-tuning of $23^ and 623^ in C7 so that the constraint from 6 ^ 57 is satisfied. 
With Mq = SOOGeV, M3 = SOOGeV, /x = 1200GeV, tan/3 = 50, and sbottom mixing S^^'^ = -0.72, the gluino 
contributions to Rg and R,[q arc shown in Fig. 2.6. R[q can reach 40% and 8% without and with the constraints 
from B XgQ and hadronic charmless B decays respectively. 
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4.3 Gluino contributions to Cq[ 

As stressed in section 2, when the CP-odd Higgs is not too heavy and tan/3 is large, neutral Higgs contribute 
significantly via Cq| ^ to the process (1). As we have stated before, only ^23^^^^^ have sizable contribution to 

C^'l 2 . The contribut ion from can be neglected due to the cancellation between contributions of the 

self energy type and penguin diagrams. Because the loop functions decease when the mass of gluino increases, 
gluino contributions to Cq| ^ in the heave gluino case are less significant than those in the not too heavy gluino 
case. So in this subsection we consider only the case where gluino is not too heavy. 

We show in Fig. 4. a the gluino contribution to Cq^ vs 5^^^ with Mq = SOOGeV, M^o = 2S0GeV, jj. = 
SOOGeV, tan/3 = 40, and sbottom left-right mixing ^fl^^ = -0.38. The other parameters are M3 = SOOGeV 
and the SU(5) gaugino mass relation at electroweak scale Mz, Mi : M2 : M3 = 1 : 2 : 7 for simplicity. Cq^ {nib) 
can be as large as 2.5, when 823^ is near 0.04. We can see that due to the constraint of Br(i? — » Xg^), 623^ is 
restricted to the order of 0.04. 

In order to analyze Cq^^, we set C'j{fni,) and Cs{mb) to 0, and keep Cg{mb) and Cio(mf,) unchanged 

from their SM values. We show in Fig. 4.6 the gluino contribution to Cq^ vs 623^ with Mq = 500GcV, 
Mao = 250GeV, M3 = SOOGeV, /x = SOOGeV, tan/3 = 40, and sbottom mixing 6^^^ = -0.38. In the figure 
the solid and dot curves (which denote the gluino contribution and all contributions respectively) are almost 
overlapping. C'q^ {nib) can be as large as 3.6, when ^23*^ is near 0.06. 

4.4 Neutralino contributions to Cfg and Cq^ ^ 

In the heavy gluino case, because we choose Mi = lOOGeV, the lightest neutralino should be light so that it 
contributes to Riq greatly and constructively with the gluino contribution, as can be seen from the dot curve 
in Fig. 1 which corresponds the case of including all contributions. We now turn to the fine-tuning case. In 
the case neutralinos can also has large effect as long as the lightest neutralino is light enough. We show in 
Fig. 3.0 the correlation of S^^-^ vs S^^-^ and Rg vs Rw- With Mq = SOOGeV, Mi = SOGeV, M2 = 300GeV, 
fj, = 1200GeV, tan/3 = SO, sbottom mixing 833^ — —0.72, Rio can reach 10%. And Rg can reach only 3%, as 
expected. We show the neutralino contributions to R'^ and R[q and the correlation 623^ vs 623^ in Fig. 3.6 
with the same values of parameters as those in Fig. 3. a. R'^q can reach the order of 6%. And i?g can reach only 
2%. 

Neutralinos can also have large contributions to Cq| ^. The neutralino contribution to Cq^ vs 623^^, with 
Mq = SOOGeV, Mao = 2S0GeV, Mi = lOOGeV, M2 ='300GeV, M3 = lOOOGeV, /x = SOOGeV, tan^ = 40, 
and sbottom mixing 633^ = —0.38, is shown in Fig. S. The Cq^^mb) can be as large as 1.5, when d^'" is 
near —0.7. It is too hard for the non-zero 823^ to generate sizable neutralino contributions to C'q^ ^, under the 
severe Br(i? ^sT) experimental constraint. 
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4.5 Wilson coefficients cfj and ^ in MSSM 

In this subsection wc switch on all the contributions of charged gauge boson , charged Higgs , chargino 
X^, neutrahno x" and gluino g and all the (5s, 62^^^, ^23^' '^23^' ^23^^' ''^^'^ ^23^^ under the experimental 
constraints. Specifically, we are interested in calculating the maximum enhancements of the Wilson coefficients 
C^, and ^ that SUSY can provide, which is important to discriminate the small tan/3 and large tan/3 
scenarios by the measurement of Br(i?s jX^ 

First, we limit ourself to the case of only one non-zero ^23^^ and including all the contributions under the 
experimental constraints. The results are presented in Figs. 1, 4 and 5. The dependences of C^g on ^'s are 
shown by dot curves in Fig. 1. It is clear that the allowed ranges of 5'^^^^^'^'^ and i?^Q are enlarged and R^(l^ 
can reach 15%. The dependences of ^ on 5'S are shown by dot curves in Fig. 4 and Fig. 5 with different 
values of gaugino masses. In the case of M3 = 500 GeV, because non-zero 5^^^ can induce scalar up-type quark 
2 — 3 left-left mixing, the chargino can have large contributions and Cq| ^ can reach about 3.5, as shown in the 
Fig. 4.«. As for the dependence on 62-^^, not like the non-zero (523"^ case, both the chargino and neutralino 
contributions are small, as shown in the Fig. 4.6. In the case of M3 = 1000 GeV, both chargino and gluino 
contributions are important, Cq^ ^ can reach about —2, and the dependence on 62^^ is quite different from that 
in the case of switching on only the neutralino contribution, as shown in Fig. 5. 

Next, we switch on all the contributions and all the Ss, §2^^, 623^, '^23^^, '^23^^, and 6^^^, under the 
experimental constraints. We perform a Monte Carlo scan of the parameter space with the ranges, |(523'^| < 1, 
|^dLL(HH)| ^ < 0.02. There are three gaugino masses Mi. 2, 3 which are free parameters in MSSM. 

For simplicity, we employ the SU(5) gaugino mass relation at the electroweak scale Mz, Mi : M2 : M3 = 1:2:7 
and take M3 = lOOOGeV in the scan. We fix Mq, ^, yl(„,d)33 with Mg = SOOGeV, fi = 500GeV, Au33 = 250GeV, 
and Ad33 = and consider two cases, small tan = 4 and large tan /? = 50. The results of the scan are presented 
in correlations, i?io vs R'iq and Cq^ vs Cq^, in Fig. 6. In the small tan/3 = 4 case, Rio can reach 12%, R'iq is 
nearly zero, and \CQ^{mb)\ and \CQ^{mb)\ are smaller than 0.1 in some parameter region due to the smallness 
of tan/3. In the large tan/? case, Rio can reach 20%, R'lo can reach 1.5%, and CQ^{mb) and CQ^{mb) can be as 
large as 3.6 in some parameter region. 

Only six Wilson coefficients, C^q , Cq^ ^, affect Br(i3s Because ^ cannot be large in the small 

tan/? case, the main contribution to Br(i3,5 ^'^fi~) comes from c|q which can have at most an enhancement 
of order of 20% compared with the Standard Model C\q. But in the large tan (3 case, if the CP-odd Higgs boson 
is not too heavy, Cq| ^ can be as large as order of onc\ which strongly enhances Br(i?s jjL^ by a factor of 
10^ — 10'^ depending on the values of parameters, as shown in next section. 



5 Numerical results for Bg — ^ ji^ ii and B — ^ Xsii^ 

In the numerical calculations, the following input parameters have been used: 

as(mz) = 0.118, as{mb) = 0.215, aemCm^) = aem{mb) = A* = -0.038, 
Kb = 0.04, Br(B ^ Xeei?) = 0.11, nif, = 0.1057GeV, mb{m^) = 2.9GeV, mfc(mfc) = 4.2GeV, 
mP"''' = 175GeV, mP°'7mP°''' = 0.29, ms, = 5.370GeV, /s, = 0.22GeV, tb, = 1.46ps. 

As can be seen from eg. (14), the branching ratio of Bg — > Xs/i+/x~ depends on the sign of C7{m,b) and, 
as mentioned before, in MSSM there exist some parameter regions, where supersymmetric contributions to 
C7 can make its sign opposite to C^^^{mb). Therefore, under the experimental constraint from the branching 
ratio of B ^ Xg'j, two separate regions for the correlation between the branching ratios of the B Xg'f and 
B — !■ Ar.5/i+/i~ are allowed. One corresponds to the case in which the sign of C7{m,b) is the same as that in SM, 
and the other corresponds to the the case when the sign of C7{mb) is opposite to that in SM, which is similar 
to the results given in ref. [15]^. Our numerical results verify the analysis. 

^In ref. [15] the neutral Higgs boson contributions to the branching ratio of Bg — + ^i"*"/*" in the large tan/3 case was missed since 
they set the mass of a muon equal to zero. 
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As an illustration, in Tab. 1 we present numerical results of Br(i?s IJ''^ IJ'~) and Br(i? Xs/U+/z~) for 
the two set of parameter values which are in the region of the parameter space where the sign of C7(mfe) is the 
same as that in SM. The other parameters which are not given in the table are M3 = lOOOGcV and the SU(5) 
gaugino mass relation at the electroweak scale Mz, Mi : M2 : M3 = 1:2:7, Mq = 500GeV, and fj, = 500GeV. 

Table 1: Wilson coefficients C^'\ g 10 Qj 3 and Branching ratios oi B ^ XgjX^ ji" and Bg IJ^'^ jJ-' ■ 

case 5^^''/5i^^/5i^^ 's^i^sik Cr/C, Cs/C, C^/C^ 

tan/3 Cio/Clo CqJC'q^ CqJC'^^ Br(B ^ X.m+A^') Br(i?, ^ m+A^") 

A 0.79/-0.046/0.0039 0.036/-0.0095 0.34/0.10 0.15/0.06 -4.4/-0.003 

4 5.2/0.002 -0.0011/0.0013 0.0008/0.0014 4.28 x 10"^ 4.29 x lO"'' 

B 0.68/-0.044/0.018 -0.037/-0.0065 0.40/0.16 0.14/0.076 -4.39/-0.00055 

50 5.5/0.036 -2.0/-2.7 2.0/-2.7 5.45 x 10"*^ 2.56 x 10"*^ 

One can sec from the table that in the case B, i. c, the large tan/3 case, Bi{Bs can be enhanced 

by a factor of 10^, compared to SM. But in the case A where tan/3 is small, Bt{Bs A*"*"/"") is the same order as 
that in SM. These two case can be discriminated at Tevatron Run II. If the observed Bv{Bs is larger 

than the Standard Model expectation value by a factor of 10 or larger and one assumes that new physics is SUSY, 
then this will unambiguously signal the large tan/3 case. As for the semilcptonic decay Br(i? Xg/i+^u"), there 
is a 50% enhancement for the values of the set of parameters in the large tan /3 case, which is closer to the central 
value of the experiment result, (2), than SM. 

6 Conclusions 

We have examined pure leptonic and semilcptonic rare B decays, Bg IJ-'^ IJ'~ and B XsiJ.'^n~, under Br(i3 — > 
Xsj) and other experimental constraints in the minimal supersymmetric Standard Model. In particular, we 
have in detail analyzed the dependence of the relevant Wilson coefficients on new flavor violation parameters, 
off-diagonal scalar quark mass matrix elements. We find that under all the relevant experimental constraints, if 
only the gluino and SM contributions are included and assuming large mixing of left-handed and right-handed 
sbottom, c['q can be enhanced by a factor of 10% in two cases, the heavy gluino and fine-tuning between 
^23'^(^2f^) and S^^^iSi^^). In particular, it is found that C['^ can be enhanced by 40% compared with SM in 
the fine-tuning case under all experimental constraints cxccipt for that from B Xgg and hadronic charmless 
B decays. When all the contributions are included and all experimental constraints are taken into account, c['q 
can be enhanced at most by 20% compared with SM. When CP-odd Higgs is not too heavy (~ 250GeV), 
and tan/3 is large 40), neutral Higgs boson penguins with gluino and down- type squark in the loop can 
significantly contribute to the bs transition and the contributions can compete with those due to the chargino 
and up- type squark loop. Comparing with the constrained MSSM, the Wilson coefficient Cio can reach a 
larger value due to the gluino and neutralino contributions, but the largest value of Cq^ ^ allowed by all the 
experimental constraints is of the same oder. From the above results, the following conclusions can be drawn: 

A. Although a 20% enhancement of the Wilson coefficient Cio, compared to SM, which is twice of that in 
CMSSM can be reached in MSSM, it alone is far from the explanation of data if Br(Ss iJb^ = 2 x 10~^ 
will be observed at Tevatron run II. 

B. Cqj 2 can reach order of one and order of 0.01 in the large and small tan /3 case respectively and consequently 

can lead to an enhancement of Br(i?s by a factor of 10^ — 10'^ in the large tan/3 case. Therefore, if 

Br [Bg l-i''^IJ'~) > 10~* is observed, there should be new physics and tan/3 must be large if new physics is the 
MSSM. 
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Appendix 

Wilson coefBcients Ci(m^), C-{myj) {i = 7, 8, 9, 10, Qi,2) of gluino contribution are as follow 

= -4xSk ^ ^^t'x- ) + )) + SG^:,- ^ (3/3(x,„,- ) + /4(X,^,- ))) 

Xi 
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+3(-G|o|o cosa + TsG^o^o sma)fcoo{xi^^o,x^o^o))] 

'^'^^^flh^flx'' ~ '^lLx'l'^l!^X°^^'^^°^^d^X°'H°X^^ 
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-0 1 mim^o 

mim^ jn^o ^la° . r,\f i \ 

+3(rpG|'o^o - G|o^o tan/3)/coo(a;j^^o,a;^o^o))] 

+ 2{Gf^* aGf^ - ^f\oGf''%)fdaoix^^^o,X^o^o,Xl-o))] 

m A. i TTi A j TTi A j Tn a. j J 3 J 



-^^^(•■pGfi - Oli, (26) 
where Vg = — f^, = — and Xjj = mf/nij, Sw = sinOw, Cw = cosOw We have check C7, 8,9,10 with ref. [15], 

"''h° '"z" 

our results agree with them except there is a minus sign difference with their box diagram. Cq^ ^ calculated 
here agree with rcf. [25]. 

The one-loop functions are normahzed to be 1, if all the arguments are equal and set to be 1. 

fi{x) = 2{l~6x + 3x'^ + 2x^ -i)x^lii{x))/{l-x)^ 

f2{x) = 2{2 + :ix -6x^ +x^ +(ixln{x))/{l - x)^ 

fsix) = 3(l-4a; + 3a;2 -2a;2ln(a;))/(2(l-x)^) 

fiix) = 3(l-a;^ + 2xln(a;))/(l-a;)3 

Uix) = 2(7-36a; + 45a;^-16a;^ + 6a;^(-3 + 2a;)ln(a;))/(9(l-a;)^) 

feix) = 2(-ll + 18a;-9a;^ + 2a;^-61n(a;))/(3(l-a;)^) 

a; In a; 
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or x^lnx y^lny z^lnz , 

tmo[x,y,z) = -31(1 - x){x - y){x -z)^ {I- y){y - x){y -z)^ {I- z){z - x){z - yY 

We follow the convension of Haber and Kane [39], to present our Feynman rules and mass matrices of 
chargino. ncutralino. squark, ami slcptons. The chargino mass matrix X is diagonalized by two matrix U and 
V, with U*XV^^ = diag(m-± , ™~±^. The neutralino mass matrix Y is diagonalized by matrix N, N*YN~^ = 
dia,g{m^o,m^o,m^o,m^o), where m^o{i = 1,2,3,4) are positive. 

The chargino (gluino, neutralino)-squark-quark couplings are as follow 

2,6 4,6 

i—l.a — 1 i—l^ct—l 
6 

+V2gs ^ -gdlT^Gf^.PL + Gf^.PR)d + H.C. 

riLb _ / ■(/■*TuL . jr* T-luR \ y^Rb _ TT TuL I/- '^b 



Gt~, = -rS, Gf^ = r^f (28) 

The chargino(neutralino)-lepton-slepton couplings are 

2,3 _ 4,6 

jC = g iii^lGll^-PLii + g ^ x"?^l(Gfj-cPL + G£'-oPh)/x 

i=l,a=l i=l,a=l 

Gfj-o = %/2(it^7V* +iiV^)r^2, Gfj_o = -\/2f^ArarJ^2 (29) 



The Z-squark-squark couplings are 



1 2 

1 /-^^ T -^^T 1 



The Z-chargino(neutralino)-chargino(neutralino) couplings are 

G^o^q = 2-^*3-^/3 ~ '^^i^-^jA^ G^o^o = ~G^o^o (31) 
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The Higgs-chargino(neutralino)-chargino(nciitralino) couplings arc 

+ig^r[G°{G^9 Pr. + G^9 Pn) + A'{G^\-Pl + G^^ Pn)]xJ 
+45x"°[G°(G^goPL + GfgoPii) + ^°(G^f-oPL + GfoloPii)]x° 

G^-^_ = —Qj^Sa + Sjj^Ca, G^_^_ = —QijSoi + SijCa 
G_--^- = Q jj^Ca -\- Sj^Sai G^-^- = QijCoi -\- SijSa 

G^f!_ = + G«f = -g,,s;3 - 5iiC;3 

Qij = \l\viiUj2, Sij = ]J^Vi2Uji 
^Lh" _ _r)"* s c!"*r rR^° — n" « ^" n 

x"x" ~ " X^x" ~ ~'^ij^a — ^ij<^a 

X°X° — " ' x'^x" ~ ^ij'^a ^ij^a 

^RG° 



x?x° ~ -WjiCfj - sp <~'^o^o - QijCfj + b^jsp 

'^^tx° ^ '3j^*^/3 ~ '^j-^C/S G^o^o = -QijSi3 + SijCff 

Qij = l{Nis{Nj2 - + Nj3{Ni2 - Natw)) 



Sij = \{N,i{Nj2-Njit^) + N,^{Na-Nat^)) (32) 
where = sin a, = cos a. The Higgs-squark-squark coupHngs are 



Gtu, = r«f (-m„,MCa^i,+^^«^r,s«)r^f/(2m^S;3)+h.c. 

Gtlu, = r»f (-m„,Mcot/35i,- + v^Al)Tf^l{2m^) - h.c. 

Gl\, = Tlf*{muilJi5ij + VuAl cot fi)Vf^l{2m^) - h.c. 

^i^d, = (-md,MCa5», - ViA%s^)V%^l{2m^cii) + h.c. 

Gf = (-md^A*Sa% + i^d^f,c„)r^f /(2m^c^) + h.c. 

^t^d, = r^f (md,/i tan/M,, - VdAf^)Tfj'^ /{2m^) - h.c. 

^£j. = rtf*{md,i^dij+VdAf.ta.nP)T$^/i2m^)-h.c. (33) 
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Figure 1: Dependence of R^q on Gluino and SM contributions are denoted by solid curves, 

and dot curves denote all the contribution. The other parameters are Mq = SOOGeV, M3 = SOOOGeV, M2 = 
1200GeV, Ml = lOOGeV, /x = 3200GeV and tan/3 = 50. 
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Figure 2: Correlations between ^j^^-^^*^^^) and and between Rg^ and R^q switching on only gluino 

and SM contributions to Rg^ and R^q in the fine-turning case. The other parameters are Mg = 500GeV, 
M3 = 500GeV, fi = 1200GeV, and tan/3 = 50. Triangle points are ruled out by S — > X^g and hadronic 
charmless B decays which require \Rs\ < 10. 
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Figure 3: Correlations between ^'^^^'^^^ ^dLR{dRL) between i?g ^ and i?Jo switching on only neutralino 
and SM contributions to i?g''' and r'iq in the fine-turning case. The other parameters are Mg = 500GeV, 
Ml = 80GeV, M2 = aOOGeV, = 1200GeV, and tan/3 = 50. 
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Figure 4: Dependence of Cg|(mfe) on S^^^^'^^^^- Gluino and SM contributions are denoted by solid curves, 

and dot curves denote all the contribution. The parameters arc Mq = 500GeV, Mj^a = 250GcV, /i = SOOGeV 
and tan/3 = 40 as well as M3 = SOOGeV and the SU(5) gaugino mass relation at the electroweak scale Mz, 
Ml : M2 : M3 = 1 : 2 : 7. 




Figure 5: Dependence of CQ^{mi,) on 623^- The solid curve denotes neutralino and SM contributions, and 
the dot curve denotes all the contributions. The other parameters are Mq = 500GeV, M^o = 250GeV, Mi = 
lOOGeV, M2 = SOOGeV, M3 = lOOOGeV, n = SOOGeV and tan/3 = 40. 
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Figure 6: Correlations between Riq vs R[q and between Cq^ vs C'q^ switching on all the contributions and all 

5|^^'s. (a) and (b) for tan/3=4 and 50 respectively. The other parameters are Mg = 500GeV, fi = 500GeV, 
Ms = lOOOGeV and the SU(5) gaugino mass relation at the electroweak scale Mz, Mi : M2 : M3 = 1 : 2 : 7, 
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